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Broad  band  vertical  and  horizontal  accelerometers  were  operated  at  the  surface 
and  at  the  1,850  foot  depth  of  a  mine  In  New  Jersey.  In  the  seismic  band  of  the 
vertical,  the  detection  threshold  on  the  surface  and  at  1,850  feet  was  Ms  ■  2.7  to 
3.2  with  a  signal  to  noise  ratio  of  2  at  25  second  periods. 

The  horizontal  was  operated  at  the  surface  and  showed  much  more  noise  than 
the  vertical.  However  the  digitally  high  passed  records  of  the  horizontal  show 
a  reduction  In  noise  of  almost  a  factor  of  10  In  the  seismic  band  (50  -  20  second 
periods).  This  allows  the  horizontal  surface  Instruments  to  have  almost  as  good 
a  filtered  noise  level  as  the  surface  vertical  accelerometers. 

A  catalog  of  all  events  taken  In  a  continuous  run  of  6  months  duration  Is 
In  progress.  All  of  this  data  Is  In  digital  form. 
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Digitally  filtered  Seismic  Instruments 
Seismic  Detection 
Wide  band  accelerometers 


UNCLASSIFIED 


LOW  FREQUENCY  DISCRIMINANTS 


During  the  course  of  this  research  contract  the  following 

developments  were  made. 

1.  A  quartz  vertical  accelerometer  was  operated  at  the  1,850  foot  level 
of  the  New  Jersey  Zinc  mine  in  Ogdensburg,  New  Jersey.  This  operation 
was  carried  out  to  compare  the  Lamont  vertical  accelerometer  with  the 
borehole  deployable  quartz  instrument.  The  comparisons  showed  a 
somewhat  lower  noise  level  from  the  quartz  instrument  in  the  seismic 
prefilter  bandwidth  (see  Technical  Report  dated  March  1,  1971).  However, 
during  the  mine  operation  a  temperature  controller  problem  was  found 

and  later  corrected  back  at  San  Diego  which  lowered  the  noise  level 
further . 

2.  The  new  vertical  quartz  accelerometer  was  operated  at  a  surfact  site 
at  San  Diego  from  January  to  July  1971.  During  this  period  data  was 
digitized  in  the  seismic  normal  mode  and  tidal  channels.  Several 
hundred  events  were  recorded  digitally.  A  student  working  for  James 
Brune  is  preparing  a  catalog  of  all  events  during  this  period. 

3.  An  analysis  of  small  teleseismic  events  was  made  to  determine  detection 
thresholds.  A  paper  on  this  is  included  (Prothero,  et  al.).  A  single 
vertical  at  San  Diego  has  a  detection  threshold  (signal  to  noise  ratio 
of  2)  for  25  second  periods  of  Hs  =  2.7  to  3.2  at  A  =  30°. 

4.  Several  large  earthquakes  yielded  earth  normal  mode  data  and  one,  the 
Colombian  event  July  31,  1970  gave  many  high  Q  overtone  modes.  The 
existence  of  the  high  Q's  of  these  modes  was  unexpected.  A  paper  on 
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this  is  included  (Dratler,  et  al.).  The  average  duty  cycle  for 
seeing  normal  mode  excitation  (20  db  signal  to  noise)  is  about  one 
day  per  month  (Block,  et  al.). 

5.  The  horizontal  accelerometer  was  placed  in  operation  in  the  same 
surface  vaul  as  the  verticals  for  several  months.  Many 
earthquakes  were  observed  and  digitized  with  these  instruments.  The 
horizontal  showed  a  much  higher  noise  level  than  the  vertical  in 
general.  However,  high  pass  digital  filtering  proved  that  the 
majority  of  this  noise  is  below  the  seismic  bandwidth  (50-20  second 
periods)  and  that  horizontal  accelerometers  can  be  operated  on  the 
surface  at  nearly  the  same  signal  to  noise  as  verticals  in  the  seismic 
band.  A  paper  on  this  is  included  (Block,  et  al.). 

6.  Preparations  were  begun  to  waterproof  the  instruments  for  use  in 
low  noise  sites  perhaps  including  mine  sites.  A  new  portable  data 
acquisition  system  was  tested. 
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Surface  Wave  Detection  with  a  Broad 
Band  Accelerometer 

We  report  preliminary  results  of  a  study  to  determine  the 
capability  of  a  broad  band  vertical  accelerometer  to  detect 
surface  waves  from  small  earthquakes.  The  instrument  uses 
a  quartz  fibre  in  torsion  and  capacitive  position  sensing1.  It 
was  designed  as  a  low  noise,  low  drift  broad  band  accelero¬ 
meter,  and  has  produced  data  on  Earth  normal  mode  excitation 
after  earthquakes  (5  to  40  cycle  tr1,  ref.  2).  For  this  study, 
the  output  was  modified  by  an  active  filter  to  copy  the  response 
of  Pomeroy  el  al.*  and  Molnar  el  al.*.  Such  a  response  does 
not  use  the  wide  band  aspects  of  this  instrument  and  the  shape 
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Fig.  I  Response  of  the  accelerometer  to  ground  motion  in  volts 
of  output  per  |im  of  ground  motion.  The  digital  fitters  shown 
are  peaked  at  20  s  and  40  s,  with  bandwidihs  of  0.02  Hz.  The 
digital  filters,  peaked  at  25  s  and  30  s,  have  almost  identical 
frequency  responses.  The  difference  in  the  widths  of  the  20  s 
ana  40  s  filters  shown  is  due  to  the  inverse  dependence  of 
frequency  and  period. 
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Fig.  2  Event  II  described  in  Table  I. 


was  dictated  primarily  by  the  shape  of  the  ground  noise 
spectrum,  which  has  a  minimum  in  amplitude  between  periods 
of  20  and  30 1.  This  minimum  may  be  inferred  from  the  known 
decrease  in  noise  amplitude  with  increasing  period  at  (tortc 
periods  and  seismic  background  noise  data*.  It  is  also  shown 
in  the  data  of  Savino  at  aL*. 

The  chief  features  of  this  ‘nstrument  of  importance  to  the 
seismic  detection  problem  (I)  a  drift  rate  which  is  small 
so  that  a  broad  band  position  detector  can  be  used;  (2)  a 
compact  design  with  internal  environmental  control  of  tern* 
perature  and  pressure;  (3)  low  ambient  and  instrumental  noise 
level,  totalling  less  than  5  x  10  ”  g*/c.p.h.  (1.8  x  10-“  cm  •-*/ 
Hz)  in  the  frequency  bend  of  interest  (/  is  the  local  acceleration 
of  gravity). 


Fig.  3  Event  S  described  in  Table  I . 


Our  data  were  obtained  at  a  curfaos  site  (vault  depth  of 
approximately  10  feet)  which  is  next  to  the  Miramar  Nara? 
Air  Station  in  San  Diego  (32.88*  N,  1 17.1*  W).  Ibis  site  is 
noisy  at  high  seismic  frequencies  because  of  jet  engine  operation 
and  traffic  on  a  nearby  freeway.  The  site  is,  however,  quiet 
at  longer  periods  and  the  instrument  response  is  sufficiently 
linear  that  no  appreciable  amount  of  high  frequency  energy 
is  coupled  down  to  low  seismic  frequencies. 

The  output  is  directly  digitized  et  1  s  intervals  and  can  be 
passed  through  narrow  band  digital  filters.  Data  taken  at  the 
Ofdensburg\  New  Jersey,  mine  site  indicated  e  noise  level 
approximately  comparable  with  that  described  by  Pomeroy 
et  al.*.  This  instrument  has  the  advantages  of  small  size,  low 
drift  rate,  and  internal  environmental  control  so  that  important 
difficulties  with  complicated  vault  construction  and  external 
environmental  control  are  avoided. 

Fig.  I  shows  the  response  of  the  instrument  to  ground  die* 
placement  at  a  given  period  in  volts  output  per  pm.  Seismic 
filter  number  4  is  the  analogue  filter.  The  unfiltered  broad 
band  response  is  also  shown  and  has  the  I  IT*  roll-off  chains 
teristic  of  an  accelerometer  with  flat  response  to  acceleration — 
T  is  the  period.  At  very  long  periods,  the  l/T*  term  become* 
small  and  the  effect  of  the  gravity  gradient  predominate*. 
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TaMe  1  The  Evam*  Described  in  Fig*.  2.  3.  4  end  S 


Event 

Epicentral  distance 

Afa  Af,  Af,  (30*)  .*  (degrees) 

h  (km) 

I.D. 

Region 

11 

3.2  3.3  4.6 

77.7 

S3 

LASA 

Tcnga  Island  Region 

3 

9 

PDE  reports  not  yet  available. 
4.8  4.1  3.9 

41.4 

43 

LASA 

Northern  Easter  Island,  Cordillera 

17 

4.3  3.3  3.2 

44.0 

33 

NOAA  PDE 
LASA 

AndreanolT  island,  Aleutians 

The  body  wave  magnitude*  (Af,)  are  as  determined  by  LASA  or  NOAA  PDE  reports  as  indicated.  The  surface  wave  magnitudes  (Aft) 
were  determined  from  our  records  using  the  surface  wave  magnitude  relation  of  Gutenberg10,  where  Afj-logio'f  + 1.66  logi0A+ 1.12.  A  is 
the  maximum  horizontal  ground  movement  (pm)  observed  for  surface  waves  with  periods  ~  20  s.  The  vertical  and  horizontal  motions  are 
assumed  equal.  Af,  (30°)  is  the  surface  wave  magnitude  of  each  event  normalized  to  a  distance  of  30’. 


Figs.  2, 3, 4  and  5  are  computer  generated  plots  of  the  seismic 
filter  output  for  several  teleseismic  events,  The  digitized  output 
has  been  digitally  band  pass  filtered  by  several  symmetric 
convolution  filters  with  bandwidths  of  0.02  Hz  and  centre 
frequencies  corresponding  to  periods  of  20,  23,  30  and  40  s. 
The  response  to  ground  displacement  for  the  20  and  40  s 
digital  filters  is  shown  in  Fig.  1.  Table  1  lists  the  magnitudes, 
depths  and  distances  of  the  events  shown.  The  surface  wave 
magnitude  of  an  earthquake  at  30*  which  would  produce  20  s 
surface  waves  of  the  same  amplitude  is  also  given  for  each 
event.  Event  number  17  of  Fig.  3  has  a  signal  to  noise  ratio 
of  approximately  3  on  the  23  s  filtered  record;  it  would 
correspond  to  an  event  with  surface  wave  magnitude  3.2  at 
30*.  The  amplitude  of  the  noise  on  the  23  s  filtered  trace  at 
this  site  (13  nr r.)  corresponds  to  the  amplitude  of  an  event  with 
a  surface  wave  magnitude  of  2.4  (20  s)  at  a  distance  of  30*. 
Other  filtering  methods  r  <ay  be  used  to  decrease  the  detection 
threshold  even  more1-*. 

An  event  with  Af,  «=  5.2  (LASA)  at  a  distance  of  77.7*  is 
shown  in  Fig.  2.  Several  body  wave  phases  are  clearly  visible; 
and  for  smaller  magnitude  events,  only  the  surface  wave 
arrivals  are  apparent.  Fig.  4  shows  an  event  with  Af,-4.8 


Pig.  4  Event  9  described  in  Table  I . 


(LASA  and  NOAA  PDE)  at  a  distance  of  41.4*  which  shows  a 
large  mantle  wave  excitation.  The  digital  filters  increase  the 
signal  to  noise  ratio  for  the  surface  wave  train  in  all  events. 
The  25  s  filter  seems  to  show  the  largest  signal  to  noise  ratio 
in  general. 

These  results  indicate  that  it  is  possible  to  achieve  a  relatively 
low  surface  wave  detection  threshold  in  a  simple  surface  vault 
in  the  presence  of  large  short  period  seismic  noise  using  the 
Block-Moore  quartz  torsion  accelerometer.  The  digital  data 
recording  system  has  also  been  found  invaluable,  as  it  provides 
a  large  dynamic  range.  Direct  digitization  also  allows  flex¬ 
ibility  and  speed  in  data  reduction.  Further  tests  will  be 
carried  out  to  determine  site  variations  in  ground  noise  ampli¬ 
tude  and  detection  threshold. 

We  thank  R.  A.  Haubrich  for  help,  and  researchers  at  the 
Ogdensburg  Mine  of  the  Lamont-Doherty  Geophysical 
Observatory,  where  the  instrument  comparisons  were  made. 
This  research  was  supported  in  part  by  the  Advanced  Research 
Projects  Agency  of  the  Department  of  Defense  and  was 
monitored  by  the  Air  Force  Office  of  Scientific  Research. 
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Fig.  9  Event  17  described  in  Table  I . 
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High-Q  Overtone  Modes  of  the  Earth 

J,  Dratlcr,  Jr.,  VV.  F„  Farrell,  B.  Block  and  F.  Gilbert 


(Received  1971  March  18  in  revised  form) 


Summary 

Overtone  modes  of  free  oscillation  of  the  Earth  have  been  observed  with 
two  instruments  at  two  separate  locations  after  the  Colombian  earthquake 
of  July  31,  1970.  Certain  overtones  persisted  with  a  signal-to-noise  ratio 
as  high  as  20  db  for  50  hr  after  the  origin  time.  Some  spheroidal  radial 
overtones  (modes  of  type  „S0)  were  found  to  have  Q's  near  1000,  about 
three  times  the  typical  Q  of  a  fundamental  mode  (0S,).  Multiplct  centre 
frequencies  and  Q's  were  measured  for  several  overtones,  with  good 
agreement  between  the  two  stations. 


Introduction 

Excitation  of  the  free  oscillations  of  the  Earth  by  large  earthquakes  has  been 
observed  for  a  decade.  The  low  ficqucncy  normal  modes  of  the  Earth  were  first 
observed  after  the  Chilean  earthquake  of  I960  (l)-(3).  Normal  mode  activity  has 
been  measured  after  several  subsequent  earthquakes  w  ith  a  variety  of  instruments  (4). 
In  these  early  investigations,  the  frequencies  and  decay  rates  of  the  fundamental 
modes  in  the  frequency  band  l-20cph  were  studied.  Spheroidal  fundamental  modes 
from  0S0  to  oSj0  were  observed  often.  Typical  Q's  for  these  fundamentals  were 
estimated  to  be  roughly  300,  except  for  0S0,  whose  Q  is  several  thousand. 

For  several  reasons,  the  spectrum  of  overtone  modes  (modes  of  type  ,S„  where 
n  >  0)  is  not  as  we>!  known.  The  overtones,  in  general,  arc  not  as  strongly  excited  as 
the  fundamentals,  and  instrument  and  earth  noise  in  the  data  hamper  their  identi¬ 
fication.  The  overtone  sequences,  having  successively  higher  and  higher  starting 
frequencies,  are  only  present  in  the  part  of  the  frequency  domain  where  the  lines  are 
dense.  Thus,  it  is  difficult  to  identify  the  overtone  modes  by  their  frequency  alone. 
Deep  fccus  earthquakes  excite  the  overtone  modes  more  strongly  than  shallow 
earthquakes,  but  large,  deep  events  are  rare. 

Alsop  A  Brunc  (5)  have  discussed  the  overtone  modes  excited  by  a  7}  magnitude 
event  at  a  depth  of  543  km.  which  occurred  on  the  Pcru-Bolivia  border  in  1963.  They 
used  phase  velocities  as  well  as  frequencies  to  identify  modes  from  a  6  hr  45  min 
record.  Second  through  fifth  overtones  in  the  frequency  range  !5-36cph  were 
observed,  but  some  of  these  were  Inter  rc-idcntified  by  Derr  (4). 

There  are  two  prerequisites  for  observations  of  free  oscillations  with  a  good  signal- 
to-noise  ratio.  One  is  a  sensitive,  low-noise  seismic  intrument  with  good  response  at 
low  frequencies,  operating  at  a  low-nois;  site;  and  the  other  is  a  large  earthquake  to 

*  Prescnl  address:  Cooperative  Institute  for  Research  in  Environmental  Science,  University  of 
Colondo,  Hould.-r,  Colorado  80 .’O’. 
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excite  the  overtone  oscillations.  In  the  present  study  two  instruments  using  extremely 
sensitive  capacitative  position  sensing,  a  modified  LaCostc  Romberg  survey  gravi¬ 
meter  (6),  and  the  Block-Moore  Quartz-fibre  accelerometer  (7),  satisfied  the  first 
prerequisite.  Fortune  supplied  the  latter,  in  the  form  of  the  Colombian  Earthquake 
of  July  31,  1970.  This  earthquake  (origin  time  1 7:08:05-4  GMT  July  31,  1970, 
epicentre  1-5°  S,  72-6°  W,  near  the  Peru-Colombian  border),  had  magnitudes 
M,  =  6-8  and  mh  =  71,  was  651  km  deep,  and  excited  the  overtone  oscillations 
strongly.  Useful  data  were  obtained  for  about  80  hr  after  the  first  arrival  with  the 
Block-Moore  instrument,  and  for  slightly  more  than  67  hr  with  the  LaCoste  instru¬ 
ment.  Overtone  modes  were  identified  from  the  decay  rates,  as  well  as  calculated 
eigenfrcqucncies,  and  good  agreement  was  found  between  the  two  instruments 
located  several  hundred  kilometers  apart. 

Data 

In  the  past  year,  numerous  recordings  of  free  oscillations  generated  by  relatively 
small  earthquakes  have  been  made  (8).  The  earthquake  discussed  here  is  the  first 
large,  deep  focus  event  we  have  observed.  Time  scries  from  the  event  were  taken 
with  the  two  different  vertical  accelerometers  at  two  locations.  The  modified  LaCostc 
gravimeter  was  located  in  Payson,  Arizona,  at  the  Tonto  Forest  Seismological 
Observatory  (34-24°  N.  Lat,  111-34°  W.  Long.).  The  Block-Moore  instrument  was 
located  at  the  IGPP  Camp  Elliott,  California,  field  station  (32-88°  N.  Lat.,  1 17-1°  W. 
Long.),  about  13  miles  north-cast  of  San  Diego.  The  distances  of  the  two  stations 
from  the  earthquake  epicentre  were  A£  *=  54-13’  for  Camp  Elliott  and  Ar  =  50-85° 
for  Payson.  Haubrich  (9)  has  studied  the  background  earth  noise  at  the  two  locations, 
using  data  from  the  same  two  instruments,  and  has  found  the  noise  level  in  the  4-20  eph 
band  to  be  about  IO'J4gJ/cph. 

Outputs  from  both  instruments  were  passed  through  active  electronic  filters  (7) 
with  gains  of  40  db  in  the  passband  from  1  to  30cph  and  40db/dccadc  rolloffs. 
Filtered  output  from  the  Block-Moore  instrument  was  recorded  digitally  on  magnetic 
tape  at  a  sampling  interval  of  I  s.  Filtered  output  from  the  LaCostc  instrument  was 
printed  digitally  on  paper  tape  at  a  sampling  interval  of  12  s,  after  additional  low- 
pass  filtering  to  avoid  aliasing. 

For  studying  the  data  from  the  Block-Moore  instrument,  a  record  beginning 
5  hr  and  ending  80  hr  after  the  earthquake  origin  time  was  used.  The  first  few  hours 
of  the  data  were  discarded  because  the  large  Rayleigh  waves  saturated  the  instrument 
electronics.  This  constituted  no  serious  loss,  as  even  more  data  was  intentionally 
ignored*  in  the  subsequent  analysis.  The  record  was  terminated  at  80  hr  after  origin 
time  as  spectra  of  records  taken  after  this  time  seemed  to  show  no  features  appreciably 
above  the  noise  level.  The  raw  data  was  interpolated  across  a  few  instrumental 
artifacts,  high-pass  filtered  to  remove  the  tides,  and  finally  low-pass  filtered  and 
decimated  to  a  Nyquist  frequency  of  60  cph. 

A  part  of  the  Fourier  power  spectrum  of  the  entire  75  hr  record  is  shown  in 
Figs  1(a)  and  2(a).  For  clarity,  the  region  from  4  to  1 2-3 cph  is  shown  in  Fig.  I,  and 
that  from  12-3  to  21  cph  in  Fig.  2.  The  theoretical  positions  of  the  fundamental 
spheroidal  modes  0S,  .arc  shown  at  the  top  of  each  figure.  With  the  exception  of 
0S2j,  all  the  fundamental  noimal  modes  in  Fig.  I  (a)  froni0S,0  to  0S2 5  can  be  identified. 
In  addition  to  the  fundamentals,  other  lilies  arc  observed  with  signal-to-noise  ratios 
greater  than  lOdb.  These  lines,  some  of  which  arc  marked  in  Figs  1(a)  and  2(a),  are 
overtone  llWilcs  of  free  oscillation.  The  most  sti  iking  example  is  the  line  on  the  low- 
frequcncy  shoulder  of  the  fundamental  0S?0.  While  it  is  well  resolved  from  the 
fundamental,  it  has  a  signal  strength  1 5  db  above  the  jitter  in  the  Fourier  spectrum, 
and  thus  is  dearly  not  noise.  From  theoretical  considerations,  it  is  identified  as  the 
tixth  overtone  of  unit  angular  order  *S,. 
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Frequency  (cycles  per  hour) 

Fig  1.  Fourier  power  spectra  from  *1  to  12  3  eph  of  data  taken  with  the  Block-Moore  vertical 
accelerometer  at  Camp  riliot,  California,  after  (lie  Colombian  earthquake  of  1970  July  31.  The 
0  dll  point  is  10* 14  g’/eplt.  (a)  Spectrum  of  tinte  se'ries  from  5  to  80  h  after  earthquake  origin  time, 
with  theoretical  positions  of  fundamental  (05,)  earth  normal  modes.  Some  overtone  lines  are 
marked  with  vertical  strokes  (b)  Spectrum  of  time  scries  from  20  to  80  h  after  origin  time,  with 
theoretical  positions(IO*  of  the  radial  oveitonej  ,S0.  Lines  arc  identified  from  theory  as  described 
in  the  text  Note  decay  of  fundamentals  relative  to  certain  ovci  tones,  (c)  Spectrum  of  time 
scries  from  30  to  80  h  after  origin  time.  Note  persistence  of  high-0  overtones. 
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Part  (b)  of  Figs  1  and  2  shows  the  Fourier  power  spectrum  of  the  record  from 
20  to  80  hr  after  the  earthquake  origin  time.  Here,  the  fundamental  modal  lines  have 
decayed  considerably,  while  the  overtone  lines  marked  in  Figs  1(a)  and  2(a)  still  have 
quite  a  high  signal-to-noise  ratio.  Only  the  lower-frcquency  fundamental  normal 
modes,  marked  in  Fig,  1(b),  are  now  identifiable.  The  theoretical  frequencies  of  the 
radial  overtones  ,S0  are  marked  in  the  figures,  while  the  observed  overtone  lines  are 
identified  as  described  below. 

In  part  (c)  of  Figs  I  and  2,  the  Fourier  power  spectrum  of  the  data  from  30  to 
80  hr  after  the  earthquake  origin  time  is  shown.  In  these  spectra,  the  fundamental 
normal  modes  have  decayed  away  almost  entirely.  The  persistance  of  certain  overtone 
lines  after  the  decay  of  the  fundamentals  indicates  that  they  have  a  Q  significantly 
larger  than  300,  the  approximate  Q  of  the  fundamentals.  In  fact,  after  m  cycles,  the 
energy  of  a  simple  harmonic  oscillator  with  quality  factor  Q  decreases  to  a  factor  of 
exp(-2nni/Q)  of  its  initial  value.  For  a  typical  fundamental  mode  with  frequency 
-lOcph  and  Q  =  300,  the  expected  decay  is  about  10  db  in  10  hr.  The  fundamentals 
shown  do,  in  fact,  have  a  decay  of  this  order,  but  some  of  the  overtones  are  considerably 
more  persistent,  suggesting  Q’s  two  or  three  times  greater  than  300.  Signal-to-noise 
ratios  for  the  final  50-hr  record  (part  (c))  are  as  high  as  20  db,  even  more  than  a  day 
after  the  origin  time. 

To  measure  the  actual  Q’s  of  the  persistent  lines,  a  4  time-lapse’  series  of  spectra 
was  made.  Segments  of  the  original  record  from  SO  to  30,  20  to  40,  30  to  50, 40  to  60, 
and  SO  to  70  hr  after  the  earthquake  origin  time  were  separately  Fourier  transformed. 
Since  the  length  of  these  sections  was  constant  and  much  greater  than  the  period  of 
any  mode  of  interest,  the  peak  of  each  line  occurred  at  the  same  Fourier  frequency  in 
each  spectrum,  and  phases  of  the  decaying  oscillations  were  unimportant.  Thus,  the 
standard  expression  for  the  decay  of  the  power  spectral  peak  from  spectrum  to 
spectrum  was  valid: 


Pk 

Pk-x 


-6>o  T 

Q 


or  lnPK 


lnP0- 


K2nf0  T 

Q 


for  Q  ►  I. 


(I) 


Here,  PK  is  the  height  of  the  Kth  ‘  time-lapse  ’  spectral  peak,  co0  *=  2 ir/0  is  the  angular 
frequency  of  the  given  mode,  and  T  =  10  hr  is  the  time  elapsed  between  the  beginnings 
of  the  consecutive  records.  The  Q  for  each  observed  high-Q  line  was  calculated  from 
a  linear  least-squares  fit  of  In  Pk  to  K  T,  according  to  the  above  formula.  The  frequency 
and  Q  measured  for  each  high-Q  overtone  line  are  shown  in  Table  I. 


Table  1 


Observed  frequencies  and  Q's  of  overtone  lines 


Mode 

Frequency  (eph) 

Q * 

Elliot 

Payson 

Elliot 

Payson 

9-0294 

90151 

672 

704 

4^0 

14-7833 

14-7819 

1173 

1156 

*5o 

17-5941 

17-5979 

938 

927 

A 

20  6690 

4$| 

7-1686 

10  3379 

10  3373 

613 

700 

•13-1541 

13-1562 

420 

380 

.0*. 

16-1592 

16-1568 

696 

574 

4^4 

8-5603 

8-5634 

4$S 

9  7276 

9-7367 

4$4 

10-8377 

10  8371 

4^7 

11-7724 

11-7806 

•  All  Q  values  have  standard  errors  between  10  and  20  per  cent  except  for  «50,  which  has  10  per 
cent  standard  error. 
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Quite  similar  results  were  obtained  with  the  data  from  the  modified  LaCosta 
gravimeter  at  Payson,  Arizona.  The  digital  record  at  Payson  began  at  0200  GMT 
August  1,  1970,  about  9  hrs  after  the  earthquake  origin  time,  and  ended  67  hr  13  min 
later.  A  Fourier  power  spectrum  of  the  data  from  10  to  70  hr  after  the  origin  time, 
after  removal  of  tides,  digital  filtering  and  decimation  to  a  Nyquist  frequency  of 
30cph,  is  shown  in  Figs  3(a)  and  4(a).  The  similarity  of  this  spectrum  to  that  of 
Figs  1(a)  and  2(a)  is  obvious.  DilTereices  in  the  apparent  excitation  levels  of  the 
various  modes  can  be  explained  on  the  basis  of  the  radiation  pattern  of  the  source. 
A  particularly  striking  example  of  this  it.  the  fundamental  node  0S2 j,  which  appears 
strongly  excited  in  the  Payson  record  and  is  non-existent  in  the  Camp  Elliott  record, 
We  infer  that  Camp  Elliott  was  nearly  on  a  node  of  the  0$22  mode  for  this  source. 

Parts  (b)  and  (c)  of  Figs  3  and  4  are  spectra  of  the  intervals  20-70  and  30-70  hr 
after  the  origin  time.  Again  there  is  good  correlation  with  the  Camp  Elliott  results. 
The  high  signal-to-noise  ratio  for  the  radial  overtone  series,  even  30  hours  after  the 
earthquake,  is  quite  remarkable.  It  can  be  ascribed  to  the  quality  of  the  site  and 
low  instrument  noise. 

A  *  time-lapse  ’  series  of  spectra  was  made  for  the  Payson  data.  Portions  of  three 
of  these  ‘  time-lapse '  spectra  are  presented  in  Fig.  3.  In  the  first  spectrum,  the  strongly 
excited  fundamentals  tend  to  obscure  the  more  weakly  excited  overtones.  However, 
as  the  series  progresses,  the  low-Q  modes  die  away,  while  the  high-Q  overtones 
persist.  The  Q’t  were  determined  for  the  lines  with  highest  signal-to-noisc  ratio,  as 
described  above.  The  results  are  given  in  Table  I .  For  those  modes  whose  Q’s  were 
reasonably  high,  there  was  quite  good  agreement  between  the  Q  values  measured  at 
Camp  Elliott  and  those  measured  at  Payson. 

The  power  spectral  density  for  all  figures  is  given  in  decibels  above  i0~>4ga/cph, 
which  is  roughly  the  Brownian  motion  noise  level  of  the  two  instruments.  It  should 
be  noted,  however,  that  the  actual  instrumental  noise  is  somewhat  higher  than  this 
theoretical  limit,  particularly  for  the  Block-Moore  accelerometer.  The  total  noise 
level  is  best  estimated  by  a  careful  examination  of  the  jitter  in  the  spectra.  For 
instance,  in  Fig.  1(a),  the  noise  level  appears  to  be  about  35  db  between  4  and  6cph 
and  less  than  30  db  above  7cph.  In  Figs  I  (b)  and  (c).  the  noise  level  estimated  from 
the  jitter  is  about  25  db.  Such  estimates  include  instrumental  noise  and  ground  noise. 

Also,  the  figures  are  not  corrected  for  the  response  of  the  electronic  and  digital 
filters  used  in  data  processing.  This  introduces  errors  in  the  absolute  sptvrtral  densities 
as  read  from  the  figures  which  are  normalised  to  the  peak  of  the  combined  filter 
response.  These  errors,  however,  are  small  and  systematic.  At  4cph,  the  figures 
read  1  db  low;  at  6cph,  the  figures  are  correct;  above  6cph,  the  figures  read  increasingly 
low  until  they  are  6db  low  at  20cph.  The  standard  error  in  the  colibration  at  the 
peak  is  ±  1  db. 

Modal  identification 

The  primary  criterion  for  the  identification  of  a  given  line  was  its  frequency.  Since 
theoretical  calculations  of  the  eigenfrequencies  have  become  both  rapid  and  accurate 
(10),  it  was  possible  to  compare  the  measured  line  frequency  to  its  calculated  value 
with  some  confidence.  When  a  given  line  occurred  at  a  frequency  well  separated 
from  those  of  neighbouring  lines,  and  when  the  observed  frequency  matched  the 
calculated  frequency  closely,  this  was  considered  sufficient  identification  of  the  line. 

Often,  however,  there  were  several  calculated  eigenfrequencies  in  the  neighbour¬ 
hood  of  an  observed  peak,  so  that  identification  by  frequency  alone  was  ambiguous. 
In  these  cases  it  was  nearly  always  possible  to  identify  the  line  unambiguously  from 
its  Q.  One  of  the  outputs  of  our  theoretical  calculation  of  the  eigenfrequencies  from 
an  earth  model  is  a  partition  of  the  potential  energy  of  a  mode  into  elastic  shear 
energy,  elastic  compressional  energy,  and  gravitational  energy.  Since  shear  is  far 
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Ft(i.  3.  Fourier  power  spectra  from  4  to  I2’5cph  of  data  taken  with  modified  LaCostc-Rombcrg 
survey  giavimete:  at  f’.tvson,  Arizona,  after  Colombian  cjithquakc  of  1970  July  31.  ThcOdd  paint 
is  10" ,J*  g-'/eph.  (a)  Speed uin  of  time  scries  from  10  to  70  h  after  earthquake  origin  time,  with 
thcoretie.il  positions  of  fundamental  (0.V,)  earth  normal  modes  Overtone  lines  are  marked  with 
vertical  stiekcs.  (b)  Speetu.m  of  time  scries  from  20  to  70  li  after  origin  ti  n;,  with  thinrctical 
positions  of  radial  o  ei  tones  „.Sa  I  ines  are  identified  from  theory  as  described  in  test.  Note  simila'ity 
to  Fig.  1(b).  (c)  Speed  uin  of  tinv:  series  from  30  to  70  h  after  origin  time.  Note  similarity  to  Fig.  1(c). 
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Flo.  4.  Continuation  of  spectra  in  Fig.  3  from  12  5  to  21  cph.  Fall-olT  with 
increasing  frequency  is  due  to  electronic  and  digital  filters  used  in  processing  the 
data.  Note  similarity  to  Fig.  2. 
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Fio.  J.  ‘Time  lapse  Fourier  power  spectra  of  data  taken  at  Pavson.  The  OdB 
point  is  10' 14  g’/eph  (a)  Spectrum  from  10  to  30  h  after  origin  time.  Note  strong 
fundamental  excitation,  to)  Spectrum  from  20  to  40  h  after  origin  time.  Note 
dominance  of  overtones,  (c)  Spectrum  from  30  to  50  h  after  origin  time.  Note  decay 
of  fundamentals  and  persistence  of  high-Q  overtones. 
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more  dissipative  than  compression  (H),  and  since  no  energy  stored  in  the  gravitational 
field  is  lost,  the  Q  of  a  given  mode  is  largely  ddennined  by  its  fractional  shear  energy. 
Modes  with  high  fractional  shear  energy  have  low  Q,  and  vice  versa.  In  cases  where 
identification  of  a  mode  was  ambiguous  on  the  grounds  of  frequency  alone,  com¬ 
parison  of  the  observed  Q  with  the  calculated  fractional  shear  energy  of  all  candidate 
modes  usually  served  to  identify  the  line. 

As  a  final  aid  in  line  identification,  a  rough  theoretical  calculation  of  the  relativo 
modal  amplitudes  was  made.  By  means  of  a  computational  shortcut  recently 
developed  (12),  it  is  possible  to  calculate  rapidly  the  spectrum  excited  by  a  given 
earthquake  source  function.  The  location  (epicentre  and  depth)  of  the  Colombian 
earthquake  was  obtained  from  the  U.S.  Coast  and  Geodetic  Survey  cards,  while  the 
fault  plane  attitude  and  direction  of  slippage  were  guessed.  The  relative  amplitudes 
calculated  from  this  source  mechanism  agreed  fairly  well  with  the  observed  relative 
amplitudes  for  the  known  lines.  These  amplitude  calculations  were  then  used  to 
corroborate  identifications  based  on  calculated  Q's  in  regions  where  frequency 
identifications  alone  were  not  sufficient. 

Results 

The  centre  frequencies  observed  at  each  station  for  the  identifiable  overtone 
multiplets  arc  listed  in  Table  I .  For  some  of  these  overtones,  the  measured  Q  is  also 
listed.  However,  this  is  done  for  only  those  lines  for  which  discrepancy  between  the 
two  measurements  is  less  than  20  per  cent. 

It  should  be  emphasized  that  not  all  the  overtone  free  oscillations  have  high  Q. 
Figs  1-4  show  many  overtone  lines  which,  although  having  greater  Q  than  the  funda¬ 
mentals,  do  not  persist  from  part  (b)  to  part  (c).  These  lines  have  higher  Q  than  the 
fundamentals,  but  less  than  the  radial  overtones.  They  have  Q's  in  the  range  500-700. 
According  to  theory,  there  are  also  many  overtones  whose  Q  is  of  the  order  of  300, 
or  the  Q  of  the  fundamentals.  These  lines  arc  not  observable,  however,  as  they  are 
initially  obscured  by  the  highly  excited  fundamentals  and  decay  at  the  same  rate,  thus 
remaining  obscured. 

A  word  of  caution  is  needed  about  the  Q  measurement.  The  representation  of  the 
normal  mode  eigenfunctions  as  spherical  harmonics  degenerate  with  respect  to 
azimuthal  degree  is  only  a  zeroth-order  approximation.  As  has  been  shown  (13), 
when  the  ellipticity  and  rotation  of  the  earth  are  considered,  the  single  lines  of  the 
zeroth  order  theory  become  multiplets  of  2/  + 1  lines,  where  /  is  the  angular  order  of 
the  spherical  harmonic.  The  existence  of  this  multiple!  structure  makes  it  impossible, 
in  general,  to  measure  the  Q  of  an  observed  mode  unless  the  individual  lines  in  the 
multiple!  can  be  resolved. 

Suppose,  for  example,  that  an  attempt  is  made  to  determine  the  Q  of  a  modal 
multiplet  by  measuring  the  line  width.  The  quantity  of  interest  is  really  the  line 
widths  of  the  individual  lines  in  the  multiplet,  which  are  related  to  the  decay  rates. 
If  these  individual  lines  cannot  be  resolved,  only  the  line  width  of  the  multiplet  as  a 
whole  can  be  measured.  This,  however,  is  not  related  directly  to  a  physical  decay 
rate,  as  it  depends  on  how  the  various  subsidiary  lines  in  the  multiplet  are  excited. 
It  is  of  interest  only  as  an  upper  bound  to  the  line  widths  of  the  subsidiary  lines.  Thus, 
in  general,  Q  values  determined  from  unresolved  multiplets  are  merely  lower  bounds 
to  the  Q  of  the  individual  lines.  A  good  example  of  this  phenomena  is  the  line  0S» 
in  Fig.  1(a),  whose  strange  shape  is  undoubtedly  due  to  differential  excitation  of  the 
47  lines  in  the  multiplet. 

Although  the  observations  of  Q  above  were  made  in  the  time  domain,  and  not  in 
the  frequency  & mr.ia,  t'ii>  reasoning  still  holds.  In  the  time  domain,  the  analogue 
of  the  phenomenon  of  line  broadening  due  to  multiplet  splitting  is  the  phenomenon 
of  beating  between  different  components  of  the  multiplet.  A  careful  consideration  of 
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the  effects  oP beating  leads  to  the  previous  result,  namely,  that  a  Q  measured  from  the 
decay  of  an  unresolved  multiplet  is  only  a  lower  bound  to  the  actual  Q  (14), 

Of  course,  the  above  considerations  do  not  apply  to  the  radial  overtones  ,S0,  as 
they  are  singlets.  As  for  the  mode*  in  the  ,S,  s  =»•»,  which  are  really  triplets,  the 
splitting  has  been  estimated  to  be  q  ■  ,e  small  for  the  modes  observed  here.  The 
effccti\e  ‘  Q ',  or /0/A/,  of  this  splitting  is  greater  t!  n  10  (4),  so  that  the  much  lower 
measured  Q  of  the  multiple!  is  a  fairly  good  estimate  of  the  individual  line  Q.  Thus, 
the  observed  Q  values  given  in  Table  1  for  the  „S0  and  „S,  nodes  should  be  rather 
accurate. 


Conclusions 

The  free  oscillations  of  the  Earth  excited  by  the  Colombian  earthquake  of  July  31, 
1970  have  \  ecn  observed  with  two  different  instruments  at  two  separate  locations.  In 
addition  to  the  fundamental  spheroidal  normal  modes,  the  data  contain  many 
spherojdal  overtunes,  some  never  previously  observed.  Analysis  of  the  data  shows  that 
the  unit  angular  order  overtones,  as  well  as  the  ,S0  radial  overtones,  have  Qt  of  the 
order  of  1000,  three  times  higher  than  the  typical  fundamental  Q.  Other  modal  lines 
of  high  Q  were  observed,  but  the  Q's  could  not  be  accurately  determined  due  to 
signal-to-noisc  problems. 

It  is  important  to  stress  that,  in  the  event  of  a  large  deep  earthquake,  these  high-Q 
free  oscillations  can  be  observed  with  any  sensitive  low-frequency  accelerometer.  In 
this  study,  the  modified  LaCoste  gravimeter  with  electrostatic  position  sensor,  a 
relatively  old  instrument,  and  the  Block-Moo  re  accelerometer,  a  new  instrument, 
were  used.  The  vital  elements  of  any  such  study  are  continuity  of  observation  and 
digital  recording  recording  of  data.  In  fact,  in  the  classic  study  of  Ness,  Harrison  & 
Slichter  (1961)  of  the  Chilean  earthquakes  of  May  I960  (1),  the  presence  and  per¬ 
sistence  of  the  high-Q  mode  4S0  was  observed,  but  it  was  not  identified.  This  was 
achieved  with  an  instrument  whose  noise  level  in  the  normal  mode  region  was  more 
than  30  db  above  that  of  our  intrsuments. 

With  the  earthquake  of  the  same  size  (M,  ss  8-2)  today,  an  enormous  wealth  of 
free  oscillation  data  could  be  obtained.  Modal  splitting  could  be  resolved  and 
compared  with  theoretical  predictions.  Many  new  overtones  could  be  observed,  and 
the  Q' s  of  these  overtones,  with  splitting  resolved,  could  be  accurately  determined. 
All  these  data,  used  as  input  to  inversion  techniques,  could  supply  a  better  picture 
of  the  internal  structure  and  dissipation  mechanisms  of  the  Earth. 


Acknowledgments 

We  would  like  to  thank  Thelma  Jones  for  her  able  and  patient  assistance  in 
typing  the  multifarious  versions  of  this  paper.  This  research  was  supported  by  the 
Advanced  Research  Projects  Agency  of  the  Department  of  Defense  and  was  moni¬ 
tored  by  the  Air  Force  Office  or  Scientific  Research  under  Contract  No.  F44620-70- 
C-0099  and  in  part  by  the  Office  of  Na\al  Research  under  Contract  No.  N000I4-69- 
A-020C-6029. 


Institute  of  Geophysics  and  Planetary  Physics, 
University  of  California, 

San  Diego,  La  Jolla, 

California,  92037 


15 


12 


J.  Dratlcr  <7  al. 


References 

(1)  Ness,  N.  F.,  Harrison,  J.  C.  fc  Slielitcr,  L.  B.,  1961.  Observations  of  the  free 

oscillations  of  the  earth,  J.  geophys.  Res.,  66,  621-629. 

(2)  Benioll,  II.,  Press,  F.  &  Smith,  S.,  1961.  F.xcitation  of  the  free  oscillations  of 

the  earth  by  earthquakes,  J.  geophys.  Res.,  66,  605-619. 

(3)  Alsop,  L.  F..,  Sutton,  G.  H.  &  Fvsing,  M.,  1961.  Measurements  of  Q  for  very 

long  period  fiec  oscillations,  J.  geophys.  Res.,  66,  291 1-2915. 

(4)  Herr,  J.  S,,  1969.  F  rec  Oscillation  observations  through  I96S,  Bull,  seism  Soc, 

Am.,  59,  2o79  2099. 

(5)  Alsop,  I  .  F.  &  Brunc,  J.  N.,  1965.  Observations  of  free  oscillations  excited  by 

a  deep  caitliquake,  J.  geophys.  Res.,  20,  6165. 

(6)  Block,  B.  &  Moore,  R.  D.,  1966.  Measurements  in  the  Farth  mode  frequency 

range  by  an  electrostatic  sensing  and  feedback  gravimetci,  J.  geophys.  Res., 
71,  4361-4375. 

(7)  Block,  B.  &  M  oorc,  R.  D.,  1970.  Tidal  to  seismic  frequency  investigations 

with  a  quartz  accelerometer  of  new  geomctiy,  J.  geophys.  Res.,  75,  1 493— 
1505. 

(8)  Block,  B.,  Dratlcr,  J.  &  Moore,  R.  D.,  1969.  Normal  modes  fiom  a  6  5 

magnitude  caitliquake,  Saturc,  226,  343-349. 

Dratlcr,  J.,  Jr.,  Normal  modes  generated  by  small  earthquakes  (part  of  Ph  D. 
thesis,  in  preparation) 

(9)  1  lattbricJi,  R.  A,,  1971.  Origin  and  ehaiaetcristies  of  microscisms  at  frequencies 

below  140  cycles  per  hour,  Monogrcphics  do  iUGGI,  in  press. 

(10)  Gilbert,  F.,  1970.  Normal  mode  cigenfrequency  calculations  (to  appear  in 

Methods  of  computational  Physics). 

(11)  Anderson,  D.  &  Archambeau,  C.,  1964.  The  auclas’  city  of  the  earth,  J. 

geophys.  Res  ,  69,  207 1  20K4. 

(12)  Gilbert,  F.,  1971.  F.xcitation  of  the  notmal  modes  of  the  Farth  by  earthquake 

sources,  Geophys.  J.  R.  astr.  Soc.,  in  press. 

(13)  Dnhlen,  F.  A.,  I96S.  The  normal  modes  of  a  rotating  elliptical  Farth,  Geophys. 

J.  R.  astr.  Soc.,  16,  329  367. 

Dahlcn,  F.  A.,  1969.  The  normal  modes  of  a  rotating  elliptical  Farth-II.  Near 
resonance  multiple!  coupling,  Geophys.  J.  R.  astr.  Soc.,  18,  397-436. 

(14)  Gilbcit,  F.  &  Backus,  G.,  1965.  The  rotational  splitting  of  the  free  oscillations 

c«f  the  Farth— 2,  Rev.  Geophys .,  3,  1-9, 


16 


REPORT  ON  A  NEW  BROAD  3AND 
VERTICAL  ACCELEROMETER 


(presented  at  ARPA  Conference 
Woods  Hole,  July  1970) 


by 

Barry  Block 


Institute  of  Geophysics  and  Planetary  Physics 
University  of  California,  San  Diego 
La  Jolla,  California 


17 


REPORT  OH  A  NEW  BROAD  BAND  VERTICAL  ACCELEROMETER 


Several  years  ago  at  the  Institute  of  Geophysics  and  Planetary 
Physics,  La  Jolla,  a  program  was  begun  to  build  a  new  generation  of 
acceleration  measuring  instruments.  This  program  has  as  its  goal  the 
creation  of  broad  band  instruments  whose  noise  and  drift  properties  are 
understood.  It  was  felt  that  the  lack  of  broad  band  accelerometers  was 
hampering  the  growth  of  some  theoretical  branches  of  geophysics  and  that 
the  narrow  band  data  gathered  was  preventing  a  clear  physical  picture  of 
many  phenomena  from  emerging.  With  due  consideration  of  the  design 
achievments  made  in  the  past,  we  felt  that  technology  had  progressed  to 
the  point  where  a  fundamental  review  was  necessary.  Our  basic  guide  lines 
were  taken  to  be 

1.  Tidal  to  seismic  frequencies  should  be  measurable  by  the  same 
instrumental  design. 

2.  Drift  and  other  noise  sources  should  be  controlled  to  a  level  where 
interesting  geophysical  information  can  be  obtained  from  tidal  to 

seismic  frequencies. 

3.  Linear  response  (i.e.,  freedom  from  non-linear  processes)  should  be 
inherent  in  the  mechanical  structure  of  the  instrument. 

4.  Internal  calibration  over  frequencies  from  DC  to  seismic  should  be  possible. 
The  advent  of  phase  sensitive  detection  allows  the  use  of  broad-band  position 
detectors  with  a  sensitivity  far  in  excess  of  what  is  needed  and  which  have 
extremely  low  detector  noise.  These  detectors  allow  a  new  design  freedom 

in  the  choice  of  the  mechanical  system.  It  is  now  no  longer  necessary  to 
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make  long  period  mass-spring  systems  to  get  the  requisite  over-ill  gain  of 
the  accelerometer.  In  fact,  there  are  many  advantages  to  making  a  simple 
mechanical  structure  with  a  free  period  from  about  1/2  to  1  second.  Such  a 
simple  structure  can  be  found  which  has  all  its  higher  modes  well  separated 
in  frequency  from  the  fundamental  mode  r.nd  which  has  small  non-linear 
coupling  from  the  higher  modes  to  the  fundamental  mode.  This  structure 
is  shown  in  Figs.  1  and  2. 

The  mass  is  cantilevered  on  a  horizontal  stretched  monolithic 
quartz  fiber  which  is  in  torsion  to  provide  the  restoring  torque.  It  should 
be  pointed  out  that  this  design  can  be  used  as  a  horizontal  accelerometer 
merely  by  turning  it  on  its  side  and  letting  the  mass  hang  down.  (Fig. 2) 

In  use  for  the  past  year  as  a  vertical  accelerometer  and  gravimeter,  one  model 
has  given  an  upper  bound  to  its  drift  rate  of  10_*°g/day.  In  other  words, 
full  scale  earth  tides  can  be  run  for  a  year  without  rezeroing.  Several 
other  models  have  shown  similar  low  drift.  This  low  drift  rate  is  achieved 
through  processes  described  in  a  recent  paper^* K 

Instrument  noise  coming  from  temperature  and  barometric  changes 
are  controlled  directly  at  the  instrument  and  not  in  a  large  vault  as 
has  been  the  custom.  Fig.  3  shows  the  mechanical  setup  of  the  instrument. 

_7 

The  mass  spring  system  is  sealed  at  10  mmHg  pressure  in  a  stainless  steel 
gold  0-ring  sealed  tank  (8  inch  diameter).  This  inner  tank  is  suspended 
from  the  lid  of  an  outer  (  9  inch  diameter)  aluminum  tank  sealed  at  1  mmHg 
by  conventional  rubber  0-rings.  In  operation,  all  adjustment  apparatus  is 
removed  from  direct  contact  with  the  inner  tank.  This  double  container 
provides  protection  from  direct  barometric  effects.  The  thermostat  is 
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wound  on  the  outer  container  and  a  thermal  feedback  system  has  been 
designed  to  provide  thermal  control  which  in  practice  does  not  drift  more 
than  a  microdegree  per  day  over  many  months  of  use. 

It  should  be  pointed  out  that  no  particular  effort  has  been 
made  to  minaturize  this  instrument  and  considerable  progress  is  possible 
if  a  smaller  version  is  needed.  The  static  calibration  and  linearization 
is  done  by  tilting  the  instrument  about  two  perpendicular  axes  and  the 
dynamic  calibration  by  using  an  electrostatic  force  applied  to  the  mass. 

The  calibration  can  be  carried  cut  to  the  percent  level  in  an  unambiguous 
way.  The  instrument  has  been  operated  simultaneously  with  three  outputs 
each  10  V  full  scale. 

_7 

a.  Tidal  channel:  Operated  at  2x10  g  full  scale.  Flat  from  DC  to 
4.8  cph,  20  db/decade  drop  after  4.8  cph. 

b.  Earth  normal  mode  channel:  Gain  of  1  from  DC  to  1  cph,  gain  of  100 

in  band  1-30  cph,  40  db/decade  drop  at  frequencies  higher  than 

_g 

30  cph;  operated  at  2x10  g  full  scale  in  pass  band  1-30  cph. 

c.  High  frequency  filter  channel:  A  variable  gain  channel  with  pass 
band  centered  at  40  second  period  whose  response  mimics  the  response 
of  the  Lamont  long  period  seismometers. 

With  this  background  in  mind,  I  would  like  to  show  you  some 

(with  a  free  period  of  1/2  second  and  Q=15) 
results  from  the  first  model/in  the  earth  mode  frequency  band  which  we 

have  somewhat  arbitrarily  called  1-30  cph.  I  will  restrict  my  remarks 

to  data  drawn  from  the  earth  mode  frequency  band  and  leave  the  seismic 

results  to  be  discussed  by  James  Brune 
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In  the  first  year  of  operation,  we  have  recorded  about  two 
dozen  magnitude  6.0  or  larger  earthquakes.  These  earthquakes  show  con¬ 
siderable  earth  mode  (1-30  cph)  excitation  which  can  extend  over  many 

.(.channel  b  above) 

hours.  Fig.  4  shows  the  earth  mode  filter  record^taken  during  the  magnitude 

(2) 

6.5,  246  km  deep  earthquake  in  the  New  Hebrides  on  October  13,  1969 

It  is  important  to  note  the  time  scale  on  this  figure.  The  background 

contains  the  earth  tides  and  we  see  that  the  earthquake  excitation  lasts 

for  nearly  twelve  hours  in  this  f”  quency  band.  A  number  of  Rayleigh  waves 

can  be  seen  and  a  detailed  look  at  an  expanded  time  scale  shows  many 

phenomena  of  interest  to  the  seismologist.  The  power  spectrum  of  the 

-24  2 

earthquake  is  shov/n  in  Figs.  5,  6  and  7.  The  0-db  point  is  1  x  10  g  /cph. 

The  ambient  noise  is  also  shown  on  these  figures  and  is  taken  directly 

from  the  preceding  record  without  any  visible  earthquakes.  This  ambient 

noise  is  the  sum  of  the  site  noise  and  instrumental  noise.  This  particular 

instrument  has  been  found  to  operate  at  its  Brownian  limit  by  cross -correlation 

with  another  instrument  but  this  cross-correlation  was  not  carried  out  at 

this  particular  setup.  The  normal  modes  begin  to  clearly  appear  at  about 

8  cph  and  show  at  least  20  db  signal  to  noise  at  frequencies  above  this. 

The  spectra  are  shown  up  to  75  cph  and  again  there  is  considerable  power 

spectral  energy  up  to  these  frequencies. 

Fig.  8  shows  a  magnitude  7.5,  20  km  deep  earthquake  at  Sumatra 
.  in  the  earth  mode  filter  (channel  b  above), 
on  November  21,  1969'  The  beginning  of  the  event  was  off  scale  at  the 

gain  used  and  the  data  collection  was  interrupted  some  18  hours  later  by 

a  site  power  failure.  There  is  considerable  excitation  still  visible 

even  after  18  hours.  The  power  spectra  are  shown  in  Figs.  9  and  10 
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(again  O-db  is  1  x  10  g  /cph).  The  signal  to  noise  ratio  exceeds  20  db 
above  8  cph.  The  power  spectra  of  many  other  earthquakes  have  been  obtained 
and  each  one  has  shown  its  own  distinctive  spectral  features, 

It  should  be  pointed  out  that  over  the  last  year,  the  duty  cycle 

for  seeing  low  order  earth  normal  mode  excitation  is  about  one  day  per 

month.  This  is  a  great  improvement  over  the  previous  rate  of  observations. 

The  large  amount  and  good  quality  of  data  gathered  will  certainly  provide 

(3) 

a  great  impetus  to  theoretical  investigation  of  source  mechanisms 

A  second  vertical  accelerometer  has  been  recently  completed  and 
is  undergoing  its  initial  tests.  It  has  several  new  features  such  as  a 
Q  =  250  for  low  Brownian  motion  noise  level  and  small  stray  capacity  to 
ground  to  increase  the  signal  to  noise  level.  This  srcond  vertical  instru¬ 
ment  is  being  operated  in  the  same  vault  as  the  first  instrument  and  a 
series  of  cross  correlation  tests  are  planned  for  the  pair. 

A  horizontal  version  of  this  basic  design  is  being  developed  by 
Jay  Dratler  and  is  contemplated  to  have  the  same  outputs  as  the  vertical 
instruments . 

In  conclusion,  the  drift,  temperature  control  and  barometric 
control  problems  seom  to  have  been  solved  for  this  instrument .  It  is 
internally  calibratable ,  linear,  and  has  a  broad-band  output.  The  first 
instrument  has  given  earth  normal  mode  excitation  for  earthquakes  of 
magnitude  6.5  and  greater  with  good  signal  to  noise.  Such  instrumentation 
can  now  provide  the  broad  band  output  necessary  for  the  next  generation  of 
seismic  investigation. 
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ILLUSTRATIONS 


Figure  1.  Exploded  view  of  torsion  spring  system. 

Figure  2.  Side  view  of  capacitor  position  sensor  and  torsion  spring  system. 

Figure  3.  Double  vacuum  can  assembly  and  tilting  platform. 

Figure  4.  Filter  output  for  New  Hebrides  earthquake  Mg  =  6.5,  h  =  246  km 
at  06:56  GMT,  October  13,  1969;  18.9°S,  169. 3°E. 

Figure  5.  Power  spectral  density  for  New  Hebrides  earthquake  of  Figure  4 
from  1-15  cph.  0  DB  on  Figure  is  1x10  24  g2/cph. 

Figure  6.  Continuation  of  Figure  5  from  12  -  40  cph. 

Figure  7.  Continuation  of  Figure  6  from  35  -  75  cph. 

Figure  8.  Filter  output  for  Sumatra  earthquake,  M#  =  7.5,  h  *  20  km 
at  02:06  GMT,  November  21,  1969;  2.1°N,  94.6®E. 

Figure  9.  Power  spectral  density  for  Sumatra  earthquake  of  Figure  8  from 
1-15  cph.  0  DB  on  Figure  is  1x10  24  g2/cph. 

Figure  10.  Continuation  of  Figure  8  from  12  -  40  cph. 
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Hg.  I  Acceleration  response  of  various  outputs  of  broadband 
accelerometer.  The  baste  broadband  response  is  that  at  the  de¬ 
tector  output.  The  other  curves  show  the  response  when  the 
detector  output  is  passed  through  three  analogue  filters.  The 
graph  shows  the  response  for  the  vertical  accelerometer.  The 
horizontal  responses  are  the  same,  except  that  the  resonance  peak 
k  slightly  shifted  and  the  relative  gain  or  the  seismic  filter  is  slightly 
different.  The  zero  dB  point  is  6.3  x  10*  V/g  for  the  vertical 
instrument,  1.95  x  10*  (V/g  or  V  rad*1)  for  the  horizontal,  and 
2.13  x  10*  (V/g  or  V  rad*1)  for  the  horizontal  seismic  filter. 


It  has  a  natural  period  of0.807  sand  a  Q  of  21.  The  horizontal 
accelerometer  is  of  essentially  similar  design.  In  this  instru 
ment,  a  simple  pendulum  24  cm  in  length,  with  a  mass  of  23  g, 
hangs  vertically  from  the  quartz  fibre.  The  fibre  stores  the 
energy  in  the  vertical  accelerometer,  but  in  the  horizontal  it 
acts  chiefly  as  a  loss-free  pivot.  Nearly  all  the  energy  is  stored 
in  the  gravitational  field,  so  that  the  natural  period  is  the 
pendulum  period.  The  results  here  were  obtained  with  an 
instrument  whose  natural  period  is  0.91  s  and  whose  Q  is  20. 

The  instruments  discussed  here  are  of  relatively  small  size 
and  feature  internal  control  of  temperature  and  pressure*.  They 


Teleseismic  Detection  with 
Wide  Band  Vertical  and 
Horizontal  Accelerometers 

Thi  observation  of  distant  seismic  events  using  a  wide-band 
vertical  accelerometer  has  recently  been  reported*.  This  is  a 
preliminary  report  of  the  observation  of  teleseisms  with  both 
the  vertical  accelerometer  and  a  new  wide-band  horizontal 
accelerometer  at  a  surface  site.  Previous  surface  measurements 
of  teleseisms  with  horizonta*  seismometers  have  been  compli¬ 
cated  by  high  levels  of  instrument  and  ground  nore1.  With 
present  instrumentation  and  the  aid  of  simple  digital  filtering, 
it  is  possible  to  obtain  signal-to-noise  ratios  for  a  horizontal 
instrument  at  the  surface  comparable  to  those  for  the  vertical 
instrument.  In  fact,  if  Love  waves  are  produced,  the  signal-to- 
noise  ratio  for  the  surface  horizontal  accelerometer  often 
exceeds  that  for  the  vertical. 

The  vertical  accelerometer,  which  has  been  described  in  detail 
already1,  uses  a  quartz  fibre  held  in  torsion  by  a  mass  of  9.5  g. 


Fig.  2  Response  to  ground  displacement  at  seismic  filter  output 
of  vertical  accelerometer  (Kale  at  left)  and  response  of  digital 
high  pass  filter  (Kale  at  right).  The  seismic  filter  output  for  the 
horizontal  instrument  is  not  readily  interpreted  in  terms  of  ground 
displacement.  The  small  bump  in  the  response  of  the  high 
pass  filter  is  the  first  sideband.  The  response  of  the  high 
pass  filter  is  cssentiall  unity  at  periods  shorter  than  45  s. 
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Fig.  3  High-resolution  cornpufr 
plots  from  digitized  records  of 
horizontal  and  vertical  seismic  Alter 
outputs.  Traces  a  and  b  show  the 
raw  outputs  of  the  seismic  filters 
for  the  horizontal  and  vertical 
accelerometers,  respectively.  Traces 
t  and  d  show  the  seismic  Alter 
outputs  from  the  horizontal  and 
vertical  instruments  after  applica¬ 
tion  of  the  digital  high  pass  Alter 
of  Fig.  2.  Although  the  scale 
is  the  same  for  all  the  traces  in 
Figs.  3-6,  the  equivalent  accelera¬ 
tion  response  of  the  horizontal 
seismic  filter  is  lower  than  that 
of  vertical  by  a  factor  of  12  (see 
•  text). 
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measure  9  inches  in  outer  diameter,  the  horizontal  instrument 
being  slightly  longer  than  the  vertical.  Internal  environmental 
control  allows  for  deployment  without  elaborate  vault  prepa¬ 
ration.  Both  instruments  are  suitable  for  deployment  in  bore 
holes. 

Both  the  vertical  and  horizontal  instruments  use  capacitative 
position  sensors  and  phase  sensitive  detection*.  Their  response  • 
to  acceleration  is  flat  from  d.c.  to  the  natural  resonance  fre¬ 
quencies  of  the  mass-fibre  systems.  As  the  natural  resonance 
frequencies  are  greater  than  I  Hz,  the  bandwidth  of  flat 
response  is  from  d.c.  to  about  I  Hz  for  both  instruments.  For 
geophysical  and  seismic  observations,  the  output  of  the  phase- 
sensitive  detector  is  passed  through  any  of  three  analogue 
filters:  a  tidal  filter,  a  normal  mode  filter,  and  a  seismic  filter. 
Fig.  I  shows  the  broadband  response  of  the  vertical  accelero¬ 
meter  at  the  detector  output  and  at  the  outputs  of  the  three 
filters.  The  curves  are  normalized  to  the  d.c.  response  at  the 
detector  output.  For  the  horizontal  accelerometer  the 
responses  are  quite  similar,  except  that  there  is  a  slight  shift  in 
the  frequency  of  the  resonance  peak  and  a  different  relative 
normalization  of  the  seismic  filter  output. 


The  dx.  gain  of  each  instrument  at  the  detector  output,  which 
is  taken  as  unity  in  Fig.  I,  has  been  carefulh/  measured.  For 
the  vertical  accelerometer,  it  has  been  determined  by  tilt 
calibration  and  by  a  fit  of  the  theoretical  Earth  tides  to  the 
tidal  filter  output.  The  result  is 

C, -(6.23  ±0.1)  x  10*  volts  lg 

where  g  is  the  local  acceleration  of  gravity.  At  this  gain, 
peak-to-peak  tidal  acceleration,  or  about  3  x  10*7/,  amounts 
to  one  tenth  of  the  full  scale  range  of  ±  10  V  at  the  detector 
output.  For  the  horizontal  instrument,  the  d«.  gain  at  the 
detector  output  has  hen  deie-nitwd  by  tilt  calibration  and 
checked  against  design  calculations.  It  is 

Gk-(l.95±0.l)x  10*  volts/g 

Although  data  from  all  thtee  analogue  filters  are  continuously 
recorded  for  both  instruments,  we  show  here  only  data  from 
the  seismic  fibers,  which  were  designed  to  match  the  response 
of  a  long-period  seismometer.  Fig.  2  shows  the  response  of 
the  vertical  accelerometer  to  ground  displacement  at  the  seismic 
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Fig.  4  Sec  caption  to  Fig.  3.  ' 
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FI*.  5  See  caption  to  Fig.  7. 


filter  output.  It  is  computed  using  the  standard  seismic  approx-  meter4  at  the  same  site.  All  the  filter  outputs  for  both  instru- 

imation  in  which  only  kinematic  acceleration  is  considered.  ments  are  recorded  digitally  on  magnetic  tape  at  a  sampling 

As  the  horizontal  instrument  measures' both  kinematic  accelera-  interval  of  I  s.  At  this  sampling  rate,  each  filter  provides 

tion  and  tilt,  interpretation  in  terms  of  ground  displacement  is  sufficient  attenuation  at  the  digitizer  Nyquist  frequency  of 

not  possible  for  it.  The  response  at  its  seismic  filter  output  is  O.S  Hz  to  avoid  aliasing  (see  Fig.  I). 

best  obtained  from  Fig.  I .  The  upper  two  traces  of  Figs.  3  to  6  are  high  resolution  com- 

Due  to  a  difference  in  the  attenuators  at  the  two  seismic  filter  puter  plots  taken  from  the  digitized  seismic  filter  data  from 

outputs.  0  dB  for  the  horizontal  seismic  filter  in  Fig.  I  is  both  instruments.  Trace  a  is  the  horizontal  seismic  filter 

2.83  *  10s  volts  per  g  of  acceleration  or  per  radian  of  tilt.  For  output,  and  trace  b  the  vertical.  The  figures  show  the  noise 

all  the  other  outputs  of  the  horizontal,  0  dB  is  1.93  x  10*  levels  at  various  times  of  night,  as  well  as  five  distant  earth- 

volts lg  as  before.  Thus,  for  the  seismic  filter  output,  the  quakes. 

equivalent  acceleration  response  of  the  horizontal  instrument  Careful  inspection  of  these  records  shows  that  the  noise 
is  a  factor  of  2.2  less  than  that  for  the  vertical.  from  the  horizontal  instrument  is  concentrated  at  frequencies 

The  two  instruments  operate  side  by  side  at  a  surface  site  slightly  below  those  of  seismic  interest.  In  fact,  most  of  the 

(vault  depth  of  10  feet)  at  the  University  of  California  Elliott  horizontal  noise  can  be  removed  from  the  records  by  digital 

Field  Station  (32.88°  N,  1 17.1 '  W).  The  horizontal  instrument  filtering.  The  lower  two  traces  in  Figs.  3  to  6  show  the  result 

is  oriented  N  43°  E  for  future  comparison  with  a  laser  strain  of  the  application  of  a  high  pass  convolution  filter  to  the  data 


FI*.  6  See  caption  to  Fig.  3. 
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of  the  upper  two  trace*.  This  high  pas*  filter,  whose  response 
is  shown  in  Fig.  2,  has  a  - 3  dB  breakpoint  at  a  period  of  4S  s 
and  is  down  20  dB  at  50  s.  After  digital  filtering,  the  pcak-to- 
peak  horizontal  noise  (trace  e)  is  reduced  by  nearly  an  order  of 
magnitude,  yet  the  seismic  signal  is  Unaffected  by  the  filtering. 
The  signal-to-noise  ratio  for  the  vertical  record  is  also  improved 
by  digital  filtering,  as  shown  in  trace  A 
Although  the  site  is  subject  to  acoustic  and  ground  noise  from 
planes  and  a  jet  engine  testing  facility  at  the  nearby  Miramar 
Naval  Air  Station,  noise  in  the  seismic  band  is  relatively  low, 
at  least  for  the  vertical  instrument.  The  horizontal  noise, 
however,  is  much  higher  and  varies  with  weather  and  the  time 
of  day.  If  the  output  of  the  horizontal  instrument  is  inter¬ 
preted  in  terms  of  acceleration,  and  if  the  difference  in  d.c.  gain 
between  the  two  instruments  is  taken  into  account,  the  vertical 
and  horizontal  noise  levels  can  be  compared.  At  night,  the 
horizontal  noise  is  a  factor  of  six  to  fifteen  greater  in  equivalent 
acceleration  than  the  vertical.  During  midday,  this  factor  is 
augmented  by  as  much  as  another  factor  of  three,  especially 
during  windy  and  variably  cloudy  days.  When  the  weather  is 
calm  and  uniformly  overcast,  as  on  foggy  days,  the  daytime 
horizontal  noise  is  as  low  as  at  night.  The  horizontal  noise  is 
due  chiefly  to  ground  tilt  induced  by  short-range  pressure  and 
temperature  fluctuations**'. 
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mental  conditions,  lead  us  to  believe  that  the  observed  noise  is 
predominantly  real  motion  of  the  Earth. 

There  are  several  features  of  interest  in  the  events  observed. 
Table  I  gives  the  identification  and  parameters  of  the  events, 
all  of  which  occurred  during  the  period  April  10*.«>  19,  1971. 
The  angle  a  is  the  smallest  angle  between  the  great  circle  path 
connecting  the  epicentre  to  the  station  and  the  line  of  maximum 
sensitivity  of  the  horizontal  instrument.  When  this  angle  is 
near  90*,  only  transversely  polarized  Love  waves  are  observed 
by  the  horizontal  instrument.  When  it  is  close  to  zero,  only 
Rayleigh  waves,  which  are  longitudinally  polarized,  are  seen. 

Event  I  had  an  angle  et  of  60*.  and  strong  Love  wave  arrivals 
were  observed.  The  signal-to-noise  ratio  on  the  digitally 
filtered  record  from  the  horizontal  instrument  was  nearly  20 
for  this  event  of  body  wave  magnitude  4.1  at  a  distance  of  10*. 
Even  after  digital  filtering,  the  signal-to-noise  ratio  on  the 
vertical  instrument  was  only  eight.  This  pattern  of  high  hori¬ 
zontal  signal-to-noise  is  repeated  in  event  7,  which  also  showed 
strong  Love  waves.  Event  7,  an  earthquake  of  body  wave 
magnitude  4.9  at  44°,  had  a  signal-to-noise  ratio  greater  than 
thirty  on  the  digitally  filtered  horizontal  record. 

Events  2,  3.  and  6  are  examples  of  earthquakes  whose 
seismic  waves  arrived  along  the  line  of  maximum  sensitivity 
of  the  horizontal  instrument.  For  these  events,  Rayleigh  waves 


Table  1  Events  Shown  in  Figs  3  to  • 


Origin  time 

Mag 

Mag 

Mag. 

M. 

Depth 

Dis¬ 

tance 

A«r 

Source 

Event 

Day  (GMT) 

(GMT)  Location 

Epicentre 

m* 

M.* 

(30*)  ♦ 

(km) 

A  (deg) 

(deg) 

1 

April  10, 1971 

10 :  43 :  40  B*ja  California 

24.3s  N,  111.0  W 

4.1 

3.0 

3.6 

33 

10 

ao 

LASAf 

2 

April  13, 1971 

05 :  17  :  38  Tonga  Island 

22.«!  S,  17S.4‘  W 

3.0 

3.0 

4.3 

33 

71 

1 

LASA 

region 

3 

April  13.  1971 

03  :  37 :  34.5  Tonga  Island 

15.9-  S,  174.0  W 

5.3 

4.5 

3.9 

73 

73 

13 

NOAA 

region 

PDET 

ft 

April  19.  1971 

02  :  43 :  52.2  Grcece-Albania 

39.0s  N,  29.3!  £ 

3.1 

4.5 

3.7 

16 

91 

13 

NOAA 

border  region 

PDE 

7 

April  19.  1971 

05 . 03 :  23  Galcpagos 

Island  region 

2.3'  S.39.6  W 

4.9 

4.6 

4.4 

17 

44 

•7 

LASA 

*  The  surface  wave  magnitude  M,  is  computed  from  the  Gutenberg  formula  <Vr,»lqt,«  -4+1.66  log,,  A+ 1  F2‘,  where  A  is  the  maximum 
zero-to-prak  amplitude  of  ground  displacement  in  pm  measured  for  seismic  waves  with  periods  near  20  s. 

t  The  surface  wave  magnitude  M,  (30*)  is  the  magnitude  of  a  hypothetical  earthquake  at  Aw  30"  which  would  produce  a  signal  of  the 
same  amplitude  as  that  of  the  given  event. 

1  Angle  a  is  the  smallest  angle  between  the  great  circle  path  from  epicentie  to  station  and  a  line  running  43*  east  of  north  (the  line  of  maxi¬ 
mum  sensitivity  of  the  horizontal  instrument). 

I  LASA*  Large  Aperture  Seismic  Array  in  Montana. 

T  NOAA  PDE*  National  Oceanographic  and  Atmospheric  Administration— Preliminary  Determination  of  Epicentres. 


Noise  originating  in  the  instrument  electronics  has  been 
measured  by  remo/ing  the  reference  signal  at  the  phase- 
sensitive  detector.  The  seismic  filter  output  for  either  instru¬ 
ment  with  no  reference  is  a  straight  l.'ne  on  the  scales  of  any 
of  Figs.  3  to  6.  Calculations  of  the  output  noise  due  to  the 
effects  of  Brownian  motion  on  the  mass  and  pendulum  have 
been  made.  The  power  spectral  density  of  fluctuations  at  the 
detector  output  due  to  Brownian  motion  is  about  3  x  I0-’ 
V*  Hz-1  for  the  vertical  instrument  and  about  I  x  I0‘*  V*  Hz*1 
for  the  horizontal.  The  horizontal  Brownian  noise  is  lower  than 
the  vertical  because  of  the  larger  mass  of  the  moving  system  and 
the  smaller  gain.  At  the  output  of  the  seismic  filter,  live  effect 
of  Brownian  motion  is  a  fluctuation  with  an  r.m.s.  amplitude 
of  18  mV  for  the  vertical  and  3.5  mV  for  the  horizontal.  As 
Figs.  3  to  6  show,  this  is  at  least  an  order  of  magnitude  down 
from  the  noise  amplitudes  observed  When  a  second  horizontal 
accelerometer  of  this  type  is  constructed,  coherence  will  be 
measured  to  determine  total  instrumental  noise  levels.  At 
present,  the  preliminary  noise  measurements  and  calculations, 
as  well  as  the  variation  of  the  horizontal  noise  with  environ- 


were  observed  on  both  instruments.  Fig.  4  of  event  2  shows 
the  characteristic  Rayleigh  wave  dispersion  on  all  four  traces. 
Under  the  assumption  that  the  horizontal  signal  derives  entirely 
from  kinematic  acceleration  and  not  from  tilt.  Fig.  4  can  be 
used  to  estimate  the  polarization  ratio  for  Rayleigh  waves. 
When  the  difference  in  gain  between  the  two  instruments  is 
taken  into  account,  the  polarization  ratio  (horizontal-to- 
vertical  equivalent  acceleration)  for  Rayleigh  waves  of  20  s 
period  is  found  to  be  about  2.6.  In  general,  the  signal-to- 
noise  ratio  for  Rayleigh  waves  is  a  factor  of  two  better  on  the 
vertical,  even  after  digital  filtering. 

Although  the  new  wide-band  horizontal  accelerometer  is 
noisy  at  this  surface  site,  the  noise  is  not  in  the  frequency  region 
of  seismic  interest.  It  can  be  removed  by  digital  filtering.  The 
digitally  filtered  horizontal  records  then  become  a  powerful 
adjunct  to  the  vertical  records  for  purposes  of  seismic  detection. 
Both  the  theoretical  power  of  Love  waves  as  a  source  dis¬ 
criminant  and  an  augmented  signal-to-noise  ratio  for  events 
producing  Love  waves  are  now  available. 

We  thank  C.  W.  Van  Sice  for  assembling  the  horizontal 
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